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Abstract 

In this paper we investigate methods to study the tt Higgs coupling. The spin and CP properties 
of a Higgs boson are analysed in a model-independent way in its associated production with a 
tt pair in high-energy e + e~ collisions. We study the prospects of establishing the CP quantum 
numbers of the Higgs boson in the CP-conserving case as well as those of determining the 
CP-mixing if CP is violated. We explore in this analysis the combined use of the total cross 
section and its energy dependence, the polarisation asymmetry of the top quark and the up- 
down asymmetry of the antitop with respect to the top-electron plane. We find that combining 
all three observables remarkably reduces the error on the determination of the CP properties 
of the Higgs Yukawa coupling. Furthermore, the top polarisation asymmetry and the ratio of 
cross sections at different collider energies are shown to be sensitive to the spin of the particle 
produced in association with the top quark pair. 



*On leave of absence from: Center for High Energy Physics, Indian Institute of Science, Bangalore 560 012, India 



1 Introduction 



One of the major goals of the LHC is to probe if the Higgs mechanism [lj-|3] is at the origin of 
electroweak symmetry breaking (EWSB). Whereas at the LHC (4j a Standard Model (SM) Higgs 
boson can be found in the whole canonical mass range and first information on its properties 
will be available, the clean environment of the International Linear Collider (ILC) jo] would be 
needed in order to determine the particle properties with sufficient precision. The interplay of the 
information available from both colliders will provide us with a clearer picture of the dynamics 
behind the creation of the masses of the fundamental particles [6] . This requires the determination 
of the spin J and the CP quantum numbers of the observed state, the measurement of its couplings 
to gauge bosons and fermions and finally the extraction of the trilinear and quartic Higgs self- 
couplings [7] to reconstruct the Higgs potential itself. 

The SM Higgs boson is predicted to be CP-even and have spin zero, i.e. J CP = ++ . Should 
extensions beyond the SM (BSM) be realised, the existence of more than one spin zero particle is 
possible. Some of these particles can have CP-odd properties or even be states with no definite CP 
quantum number. Supersymmetric theories |8, 9j for example require the introduction of at least 
two complex Higgs doublets. In the Minimal Supersymmetric Extension of the Standard Model 
(MSSM) [3 ,10 this leads after electroweak symmetry breaking to five physical Higgs states, out of 
which three are neutral, two CP-even and one CP-odd. Should CP violation beyond the one of the 
SM be provided by the Higgs sector, this would imply Higgs states with mixed CP properties. CP 
violation is one of the conditions to explain the Baryon Asymmetry in the Universe 111]. A study 



of the CP properties of the Higgs boson can hence already give some clues to BSM physics 12 



Such a study at the current and future colliders fl3] would have to establish the CP properties of 
the observed spin particle(s) and determine the amount of CP-mixing in case of CP violation. 
This can be done by establishing the tensor structure of the Higgs couplings to Z or W pairs and of 
the ff& couplings, where <3? denotes a general J = state. Any deviation from the SM prediction 
can then be interpreted in the framework of a given model. 

In the SM, the largest Higgs boson couplings are the ones to the heaviest fermions and to a 
massive W or Z gauge boson pair. Therefore Higgs production via these couplings or its decays into 
a heavy // or WW, ZZ pair provide the best processes to obtain information on its spin and parity. 
The couplings to a pair of photons or gluons, which proceed through loops of these particles, can 
be used, too. A plethora of observables, constructed out of the kinematic variables of the different 
particles involved in the production and decay of the Higgs boson, can provide this information. 
The strategy for a model-independent test and analysis is to assume the most general form of the 
relevant Higgs couplings consistent with Lorentz invariance and gauge invariance and explore how 
these couplings can be constrained by using the mentioned observables. This is usually possible by 
exploiting properties of kinematical distributions such as threshold effects and angular correlations 



in Higgs decays into W or Z pairs or production in 77 fusion, without 12, 14-18 and with the 



inclusion of CP-violating effects [19j-[27j . In the decay channel into a Z boson pair even small signal 
samples, as available at the moment of Higgs discovery at the LHC, might be sufficient to draw 



first conclusions on the CP properties of the Higgs boson 27 . The Higgs production processes 



through Higgs-radiation and gauge boson fusion in e + e collisions also provide observables to 
extract spin parity informations of the Higgs boson as well as the couplings of the Higgs boson with 



a gauge boson pair 16,28-30 , with CP violation included in [19 31-331. In addition, the azimuthal 



angular distribution of the two outgoing forward tagging jets in weak boson fusion and in gluon 
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fusion initiated processes can be used 34 -36 . In vector boson fusion the CP nature of the Higgs 
VV (V = W, Z) coupling is probed, whereas gluon fusion provides this information for the Higgs 
coupling to a top quark pair. It is found that a CP-even Higgs boson can be clearly distinguished 
from a CP-odd Higgs state. 

The coupling of a pseudoscalar Higgs boson to a VV pair is always loop induced and hence 
suppressed with respect to the tree level coupling of a scalar state. As a result, the best among 
the different methods suggested to study the CP-mixing are those which use the couplings of a 
J = state to a pair of photons 14,21-25] at the photon collider option [37] of the ILC, and to 



the heavy fermions t or t |38f|47| at the LHC and the ILC. The advantage is that these particles 
couple democratically to the CP-even and the CP-odd components of the Higgs boson. As the 
photon collider option is the most remote one to be realized, the study of the associated production 
of a J = state with a // pair at the LHC or ILC or the analysis of the decays of a heavy Higgs 
produced via gluon fusion into a // pair are the most promising candidates to get unambiguous 
information on the CP properties of a state with indefinite CP quantum numbers. 

The top quark is the heaviest fermion discovered so far and, providing the largest Yukawa 
coupling, plays a special role. At a future ILC Higgs boson production in association with a top 



quark pair, e + e — > tt& 48-52 leads to sufficiently high rates 53 and can hence be used to 
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extract CP information 41,43-45 by exploiting angular correlations and/or the polarisation of 



heavy fermions [44 ,45 . These can also be exploited in the Higgs decays to heavy fermions 38^ 
Reference [43] employs an optimal observable technique to show that the discrimination between a 
CP-even and a CP-odd case at a high level of statistical significance is possible with low luminosity, 
and that a non-zero CP-violating coupling can be determined at the la level with J C = 100 fb _1 . 
In Ref. [44| the extraction of the Higgs CP quantum numbers in ii<3? production at a future ILC 
through the measurement of the total cross section and its energy dependence as well as of the 
top polarisation asymmetry has been discussed. In the current paper we extend the analysis in 
two ways. Firstly, we allow for more than two neutral spin states while writing the most general 
CP-violating //$ vertex. Secondly, we investigate possible sensitivities to probe this general //<& 
vertex in the combined study of three different observables. The observables we use are the total 
cross section a and its energy dependence, which has been shown to exhibit an unambiguously 
different behaviour for CP-even and CP-odd Higgs bosons, the top quark polarisation asymmetry 
Pt and the up-down asymmetry of the antitop quark with respect to the top-electron plane. 
The latter can directly probe CP violation. We also discuss how well the energy dependence of 
the cross section and the top polarisation asymmetry can differentiate between the case of a spin 
1 particle produced in association with a ti pair and that of a spin particle. 

The outline of the paper is as follows. In section 2 we will discuss the observables and their 
individual sensitivities to the Higgs CP properties. Then the combination of all observables will be 
investigated with respect to an improvement in the determination of the CP quantum numbers. In 
section 3 the associated production of a spin 1 particle with a top quark pair will be discussed to 
investigate to what extent the total cross section and the top polarisation asymmetry can help to 
distinguish a spin 1 from a spin state. We will conclude in section 4. 
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2 The observables 



In this section we will discuss the observables for the production of an arbitrary CP- violating Higgs 
boson $ in association with a top-quark pair at a future e + e~ linear collider, 

e+e~ -> . (1) 

The value of the total cross section a and its energy dependence exhibit sensitivity to the CP nature 
of the Higgs boson and allow for a distinction between a CP-odd and a CP-even Higgs boson. The 
top quark polarisation asymmetry Pt serves as a further observable. Both a and Pt are CP-even 
and hence do not test CP violation. The tiQ production process, however, exhibits already at 
tree-level CP violation in the interference between the diagrams where $ is emitted from the top 



(antitop) and where it is emitted from the Z boson 41 . The up-down asymmetry A$ projects out 
this interference term and thus tests CP violation. In the following we will discuss the observables 
a, Pt and in detail as well as their sensitivity to the Higgs CP properties. We will examine if 
the longitudinal polarisation of the initial beams can help to improve the sensitivity. For the 
degree of polarisation, the standard ILC values of P e - = —0.8 and P e + = 0.6 will be assumed. 
Positive values of P e - , P e + correspond to right-handed polarisation. Finally, we will combine all 
observables to derive the errors on the CP parameters in a general CP-violating ii<3? coupling by 
performing a x 2 test. 



2.1 The total cross section 

The associated production of a SM Higgs boson with a top quark pair, e + e~ — > ttH, can be 
measured with an accuracy of 0(10)% for Mh ^ 200 GeV [53] . At tree level [48], it proceeds 
through the Higgs radiation off the t, i lines and a diagram with the Higgs boson produced in 
association with a Z boson which then splits into a tt pair, cf. Fig. [T] This diagram only 
contributes with a few percent as long as y/s < 1 TeV. In fact, the bulk of the cross section stems 
from the splitting of the virtual photon into ti. 




Figure 1: Feynman diagrams contributing to the associated production of a Higgs boson $ with a 
top quark pair. 



In models with more than one Higgs boson as e.g. a general 2-Higgs doublet model (2HDM) 
or supersymmetric extensions of the SM, there are additional diagrams with a CP-odd (even) 
Higgs boson splitting into ti for (pseudo-)scalar Higgs production. In CP-violating models both 
Higgs bosons are CP-mixed states, one of them splitting into a top quark pair. This contribution 
is in general small unless Higgs decays into ti are kinematically allowed. We assume here that 
by the time the CP quantum numbers of Higgs bosons will be tested, all Higgs bosons will have 
been discovered and their masses will have been determined. By applying appropriate cuts on the 
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invariant mass of the tt pair we can hence safely neglect these additional diagrams^] 

In order to discuss a general CP-mixed Higgs state we parametrize the tt<& coupling by 

C m = -i . a ^~( a + = -WttHia + i&7s) , (2) 
sm Ow 2,M\v 

where 9\y denotes the Weinberg angle, m t the top quark mass and My/ the W boson mass. The 
coefficients a and b are assumed to be real and in the SM are given by a = 1, b = 0. A pure 
pseudoscalar coupling is provided by a = 0, 6 / 0. A coupling to a Higgs state with indefinite 
CP quantum numbers is realized if simultaneously a / and 6/0. The exact values of these 
parameters depend on the model under consideration. In the MSSM for example the Higgs couplings 
to up-(down-)type quarks are suppressed (enhanced) for large values of tan/3, the ratio of the 
vacuum expectation values of the two complex Higgs doublets [3 54 . From the experimental side 



the upper bounds on the electric dipole moments of the neutron and the electron provide important 
constraints on non-standard Higgs sector CP violation [12]. Within a general 2HDM with maximal 



CP violation, values with \ab\ < 2 are in accordance with low-energy constraints [40,55 . We work 
in a model-independent approach and in the following let a and b vary between —1 and 1, if not 
stated otherwise. 

For the calculation of the total cross section the ZZ& coupling will be needed, too. It is 
parametrized in terms of the SM coupling gzzH by a parameter c, i.e. 

The ZZ<& coupling and hence the parameter c will be determined from other channels [4] [5] [56] [57] 
so that the analysis for the determination of the Higgs CP properties can be performed for a fixed 
value of c. In order to reduce the number of free parameters, we will choose c = —a in our analysis. 
This is a reasonable choice as c = in case of a CP-odd Higgs boson and |c| < 1 in generaQ For 
a = 1 the ZZQ coupling then takes the SM form. In principle, ZZ can also couple to the CP-odd 



part of the Higgs boson [5&-61 . This coupling is, however, zero at tree- level and only generated 
through tiny loop corrections. Therefore we consider only ZZQ couplings without contributions 
involving the CP-odd part of <£. 

Neglecting the small contribution of the diagram involving the ZZQ vertex, the differential 
cross section with respect to the scaled energies x\ 2 = 2-E t f/y / s, can be cast into the form 



da _ a 2 J 

dx\dX2 47TS 1 



QeQt+ {i-KY + l-h z 



(a 2 Ff + b 2 Ff) 



where a~ l =a~ 1 (s) ~ 128, h z = M\js and Vf = (2I^ L - AQ jsly) / (Aswcw),CLf = 2lj L / (Aswcw) are 
the usual Zff couplings given in terms of the fermion charge Q f and the third component of the 



In case the Higgs bosons are not discovered because of suppressed couplings to gauge bosons or too heavy masses, 
the process is not affected by these additional Higgs bosons either. Their influence is negligible in the first case, and 
they won't play any role in the second case at sufficiently low cm. en ergy. 

2 This follows from a sum rule for the Higgs gauge coupling |58[|59| in models with a Higgs sector containing only 
Higgs singlets and doublets. 
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Figure 2: Total cross section in fb for SM (full) and purely CP-odd (dashed) Higgs production in 
association with a ti pair as function of the cm. energy for M$ = 120 GeV and unpolarised 
beams. The pink (grey) band indicates the statistical error at ±5<7 with J £ = 500 fb -1 . 

weak isospin (f = t,e). We have introduced the short-hand notation sw = sin 8\y, cw = cos 9\y- 



The expressions for the form factors for scalar and pseudoscalar Higgs bosons are given in 



Refs. 48 . Note that in the spin averaged matrix element terms proportional to the combination 
ab of the CP-even and CP-odd coupling parameters in the Yukawa coupling do not appear. 

The total cross section is a CP-even observable and not sensitive to possible CP violation. Its 
threshold rise, however, is strikingly different for the scalar and pseudoscalar case as has been 
shown in Ref. [44] . Parametrizing the deviation p from the threshold by 

P= l 7T-^r> ( 5 ) 

V s V s 

the latter shows a softer dependence on p, which is given by ~ p 3 . The rise for a purely CP-even 
Higgs boson H on the other hand is ~ p 2 . The threshold rise can hence be exploited to distinguish 
a pseudoscalar from a scalar Higgs boson. Taking into account only statistical fluctuations, for a 
120 GeV Higgs boson an integrated luminosity of 500 fb -1 is sufficient to distinguish the purely 
pseudoscalar from the SM case at 5a confidence level, as can be inferred from Fig. [2} According 
to our coupling parametrization, here and in the following SM refers to a = 1,6 = in the tt$ 
coupling and purely pseudoscalar refers to a = 0, b = 1, if not stated otherwise. 

In the following we discuss to what extent the value of the total cross section can be exploited to 
extract the CP properties of the Higgs boson. The cross section is subject to higher-order (SUSY-) 



QCD 49 -51 and electroweak (EW) [52] corrections. The QCD corrections can be significant for 
collider energies near the threshold. In the continuum, for ^/s = 1 TeV they are of moderate size 
for both the scalar and pseudoscalar Higgs boson. The EW corrections can reach about 10%. For 
a precise determination of the top Yukawa coupling the corrections have to be taken into account. 
In this analysis we neglect their influence in a first approximation. We will supplement our analysis 
later with the use of the polarisation asymmetry of the top quark which we expect to be less 
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Figure 3: Contour plot of the total cross section cr(e + e~ —> tt$) in fb in the a — b parameter plane 
for Mh = 120 GeV and \fs = 800 GeV with unpolarised (left) and polarised (right) beams. 
The grey code indicates the value of the cross section. 

sensitive to such corrections. 

The total ti§ cross section values are shown for unpolarised beams in Fig. [5] (left) as contours 
in the a — b plane. We have chosen M<j> = 120 GeV and the cm. energy y/s = 800 GeV. Note 
that here and in the following the diagram with the ZZ& coupling is always taken into account 
in the numerical analyses. Since the cross section depends quadratically on a, b (for our choice of 
c = —a) the contour lines are symmetric with respect to the sign of a and b so that only results 
for positive values of a and b are shown. The cross section decreases with decreasing values of a, b. 
The SM value of the cross section amounts to 2.78 fb. The cross section for a purely CP-odd Higgs 
boson is with 0.4 fb much smaller than the one for a SM Higgs boson, see also Fig. [2] Furthermore, 
the dependence on b is rather flat in contrast to the dependence on a. We hence expect the cross 
section to be mostly sensitive to a. Polarisation of the initial e^- beams can help to increase the 
cross section for our choice of P e ± , as shown in Fig. [3] (right). It increases the cross section by about 
a factor of 2, independent of the CP nature of the Higgs boson. In the Appendix the replacements 
are listed, which have to be applied in the formula of the cross section in case of initial beam 
polarisation. 

Figure [4] shows the errors on a, b extracted from the total cross section, and hence the sensitivity 
of this observable to the Higgs CP parameters. For the error determination we define a level of 
confidence / to identify the area in the a — b plane for which the value of the observable 0(a,b) 
(here a(a, b)) cannot be distinguished from the reference value 0{oq, bo) at a specific point (ao, bo). 
Ignoring systematic errors, it is given by 



where AO(ao, bo) is the statistical fluctuation in O at an integrated luminosity £ chosen to be 500 
fb -1 in the following, if not stated otherwise. For the cross section it is given by 



\0(a,b)-0(ao,b )\ = fAO(a ,b ) 



(6) 




(7) 
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Parameter a Parameter a 

Figure 4: Errors Aa + (upper left) and Aa~ (upper right) on a as well as A6 + (lower left) and 
Ab~ (lower right) on b, extracted from the total cross section a(e + e~ —> ti<&) at la confidence level 
for Mj, = 120 GeV at ^/s = 800 GeV with / C = 500 fb -1 . The e ± beams are unpolarised. The 
colour code indicates the magnitude of the respective error. 

The errors on a and b for the point (ao,&o) are then determined by the maximal extensions Aa + 
(Aa~) in positive (negative) a direction and Ab + (Ab~) in positive (negative) b direction, necessary 
to reach the area outside the thus defined range of insensitivity. 

As anticipated from the strong dependence of the cross section on a and as can be inferred 
from the figure, the cross section is very sensitive to a. In most of the parameter range a can be 
determined with an accuracy of 0.2 or better. The error Aa~ is large for a < 0.4. This can be 
understood by realizing that the area of insensitivity at la for a point (ao,&o) i n this parameter 
region is given by an elliptic band in the (a, b) plane around this point. Since the value of b is 
not known a priori from any other measurement, the area outside the band of insensitivity is only 
reached for large Aa~. An illustrative example is shown in Fig. [5] (left). For a > 0.4 the insensitive 
area cuts the edge of the parameter space \b\ < 1, so that the error Aa~ is smaller, cf. Fig. [5] 
(right). The sensitivity to b is only good where the influence of the dominant contribution to the 
cross section from terms proportional to a is small. This is obviously the case for small a values. 
The total cross section can therefore be exploited to determine the CP-even part of a ti& coupling. 
From the discussion concerning Aa~ we expect that the sensitivity on a will significantly improve 
if b can be constrained from other observables. We furthermore found that polarised beams 
improve the sensitivity only marginally. 
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Figure 5: la band of insensitivity for (ao,&o) = (0.2,0.3) (left) and (ao,&o) = (0.7,0.3) (right). 
The corresponding errors Aa + , Aa" are indicated in the figures. The errors Ab^ are derived 
analogously 



2.2 The polarisation asymmetry of the top quark 



Since the top decay width is large (I\ ~ 1.5 GeV) the top quark decays much before hadronization. 
Its spin information is hence translated to the distributions of the decay products and not contam- 
inated by strong interaction effects. As the lepton angular distribution in the decay t — > bW —> blv 
is not affected by any non-standard effects in the decay vertex, it is a pure probe of the physics 
associated with the top quark production process 22,62 . Thus the polarisation asymmetry Pt of 



the top quark is another observable expected to probe the Higgs CP properties. Both for polarised 
and unpolarised beams it is given by 



a(t L ) - o(t R ) 



(8) 



where t^R denotes a left-, right-handed top quark. Note that, since Pt is a CP-even quantity, 
the polarisation asymmetry of the anti-top is the same as Pt but with opposite sign. The cross 
section for an L-,i?-polarised top quark can be decomposed into a term oo proportional to the 
spin-averaged cross section and a term a\ proportional to the helicity of the top quark, 



&l,r = Co — j A = — 1, +1 for L, R , 
with the total cross section given by 

cr(W$) = 2cr . 

We can hence express Pt as 

2(71 



Pt 



2a ' 



(9) 
(10) 

(11) 



where the differential cross section corresponding to a±, neglecting the small contribution of the 
diagram involving the ZZ& vertex, is given by 



dai a 

dx\dX2 47TS 



a t9ttH 



QeQtVe (aj + vj)vt 

i-h z {\-h z y 



i 2 G H + b 2 G A ) 



(12) 
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with the form factors 

G H (xi,x 2 ) = P(2 — xi — x 2 )$xi(l — xi)(xi — x 2 )(l — x 2 ) — 16/if (2 — x\ — x 2 ) 

—2ht[4 + x\(-3 + xi)(2 + xi) — 6x2 + 8x\x 2 + (1 — + hq>[x\{x2 — x\) 

+Ah t {2 - Xl - x 2 )]) (13) 

and 

G A (xi,x 2 ) = P{x\ - x 2 ){2 - x\ - X2)|/i$(4/it + xi) + (1 - xi)[xi(l - X2) 

-2h t (2 - Xl - x 2 )]} . (14) 
We have introduced the prefactor 

P = . 1 (15) 

v /l-^(l-Xi)2xi(l-X 2 )2 

and h$ = M^/s, ht = ruffs. Since Pt is given by a ratio of cross sections, this quantity has the 
advantage of being independent of an overall model-dependent normalization of the ti§ coupling. 

The total cross section a and the helicity projected cross section 2o\ are shown in Fig. [6] (left) 
and (right), respectively, as a function of ^fs for a scalar SM Higgs and a purely pseudoscalar 
Higgs and for two different Higgs masses M$ = 120 and 150 GeV. Like in the case of the total cross 
section, a\ shows a flat threshold rise for the pseudoscalar and a steep one for the scalar Higgs. The 
dependence on p, which parametrizes the small deviation from the threshold cm. energy, differs by 
approximately one power for the CP-even and the CP-odd case, like in the case of a = 2<tq. This 
results in a threshold rise for Pt which is approximately the same for the scalar and the pseudoscalar 
Higgs as can be inferred from Fig. [7j Due to the reduced phase space the cross section and the top 




600 800 1000 1200 1400 1600 1800 2000 600 800 1000 1200 1400 1600 1800 2000 

Vs [GeV] Vs [GeV] 

Figure 6: The total cross section a (left) and the helicity projected cross section 2<7i (right) for a 
SM and a purely pseudoscalar Higgs boson as a function of -^/s for two masses M$ = 120 and 150 
GeV with unpolarised beams. 
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Figure 7: The top polarisation asymmetry for a SM and a purely pseudoscalar Higgs boson as a 
function of yfs for two masses M$ = 120 and 150 GeV with unpolarised beams. 



polarisation asymmetry decrease with increasing Higgs mass. 

Furthermore, the values of ctq, g\ and of Pt are very different depending on the CP quantum 
numbers of the spin state. The difference is solely due to different form factors F^ 2 , F^ 2 (G H , G A ) 
for the scalar and pseudoscalar Higgs in gq {g\ ) . There is no difference in the coupling structure for 
the scalar H and the pseudoscalar A, cf. Eqs. Q and (12). At very high energies, the chiral limit 
is reached and the form factors for H and A (and hence the corresponding values of o"o and g\) 
become equal for a scalar and a pseudoscalar Higgs, up to the contributions from the diagram with 
the ZZ§ coupling, which are subleading at high cm. energy. This behaviour can be anticipated 
from Figs. [6] (left) and (right). For the polarisation asymmetry, being the ratio of gq and G\, the 
chiral limit is reached much more slowly. We also studied this behaviour of Pt at high energies in 
the process e + e~ —> &6<3? which, apart from the couplings and masses, is identical to ti& production. 
Due to the smaller final state masses, however, the process takes place at lower energies and the 
chiral limit can be observed at much lower energies compared to the ti<& final state. 

The values of Pt are shown as contour plots in the a — b parameter plane in Fig. [8] for M$ = 120 
GeV, y/s = 800 GeV and unpolarised beams. The top polarisation asymmetry in case of a 
purely pseudoscalar state is with 0.03 much smaller than the one for a SM Higgs boson with 0.11. 
Furthermore, Pt rapidly increases with a for large values of 6 and small values of a and then 
approaches a nearly constant behaviour in a. For small 6, Pt is constant in a (due to our choice 
c = —a), as the 6 contribution to gq and G\ does not play a role and can be neglected. The coupling 
factor a 2 hence cancels in the ratio of g\ and gq so that Pt is simply given by the corresponding 
ratio of the form factors and is independent of a. The same of course holds for small values of a. 
Here Pt is almost constant in 6. Overall, the change of Pt with 6 is small. The largest values of Pt 
are reached for large a. Polarising the beams increases Pt independently of the CP nature of 
the Higgs boson by roughly a factor of 3. The replacements to be made in the formulae of gq, g\ 
for polarised beams are given in the Appendix. 

The errors on the extraction of a and 6 from the top polarisation asymmetry have been calculated 
as well, with AO(ao,6o) of Eq. ^ given by 

1 



AO(a , b ) = AP t = ~^=Jl - P? ■ (16) 
VgC v 
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Figure 8: Contour plot of the top polarisation asymmetry Pt in the a — b parameter plane for 
M H = 120 GeV and = 800 GeV. The e ± beams are unpolarised. The grey code indicates the 
value of Pt. 



The errors on a and b for unpolarised beams turn out to be larger than 1. Polarisation only slightly 
improves the situation. Whereas the observable Pt can be used to distinguish a pseudoscalar from 
a scalar Higgs, the use of Pt alone does not allow for an accurate extraction of a and b, unless the 
parameter range has been constrained elsewhere before. 

2.3 Ratios of cross sections 

Scalar and pseudoscalar Higgs production differ by their threshold behaviour. A simple method 
to test this different behaviour is the measurement of the total cross section at two different cm. 
energies. These should both be below ~ 800 GeV, so that they lie within the region of rising a, 
and the number of events should not be too small. As cm. energies we therefore chose ^/s\ = 800 
GeV and ^/s2 = 600 GeV and investigated the ratio 



The advantage of investigating ratios of cross sections is that the effect an overall model-dependent 
normalization of the top quark Yukawa coupling and some systematic errors in the measurement 
can be avoided. We found, however, that the results for Aa^ and Ab^ show the same behaviour 
as for the observable Pt. The ratio R is a good observable to distinguish a purely CP-odd from 
a CP-even Higgs boson in case of a CP-conserving Higgs sector, or to distinguish a CP-mixed 
Higgs state with a large CP-odd contribution from a CP-even Higgs. It is difficult, however, to 
extract the absolute values of a and b in a model-independent approach over the whole considered 
a — b parameter range. In the following we will therefore not use R, but instead Pt (together with 
the total cross section and the up-down asymmetry) to determine the CP quantum numbers. It 



R = 




(17) 
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may be possible to measure the ratio R more accurately, but this requires running at two different 
collider energies, whereas for the measurement of Pt, while more complex, running at a single 
energy suffices. 



2.4 The CP-violating up-down asymmetry 

Contrary to the total cross section and the top polarisation asymmetry the up-down asymmetry 
Aff, of the antitop with respect to the top-electron plane is an observable which is sensitive to CP 
violation. Defining by q a ^ the four- momenta of the incoming e~ , e + and by the four-momenta 
of the top, antitop, respectively, the angle <fi between the antitop direction and the top-electron 
plane is given by 

sin0= m^m (18) 

\P2\\qa X Pl\ 

with the totally antisymmetric Levi-Civita tensor e. The up-down asymmetry of the ti& cross 
section a is defined as 

o-(up) - a(down) 

A <t> = ~i — r~i — Ta — V ' V 19 J 

cr(upj + a (down) 

where 'up' ('down') denotes the value of the cross section if the integration over (j) is performed for 
4> € [0,7r) {(f) € [tt,2tt)). It turns out that is given by the interference of the diagram where 
the Higgs is radiated from the top-quark with the diagram where the Higgs is radiated from the Z 
boson [41| . The asymmetry can be cast into the form 

A, = ^ , (20) 
a 



with the numerator of Eq. (20) denoting the asymmetric interference term proportional to 6c, 



which can be derived from the differential cross section given in Eq. (29) of the Appendix. The 
denominator is given by the total cross section a = <r(up) + cr(down). A non- vanishing A$ is an 
unambiguous indicator of CP violation. As it is proportional to b, the ti<& coupling must contain 
a CP-odd part. If, however, the Higgs sector is CP-conserving and a purely CP-odd Higgs boson 
is radiated from the top line, the coupling to the Z boson must vanish at tree- level, i.e. c = 0, so 
that in the CP-conserving case A^ = 0. The contour lines of A^ are shown in Fig. [9] in the a — b 
parameter plane for M<j, = 120 GeV and yfs = 800 GeV. The incoming e^ beams are unpolarised. 
As we put c = —a, the up-down asymmetry depends in the numerator linearly on o and b. In the 
denominator it depends via the total cross section quadratically on a and b. The asymmetry is 
shown for positive values a and b. It can be inferred for the negative values from its asymmetric 
behaviour with respect to a sign change of a or b. We found that the polarisation of the initial 
beams does not help to increase Aj,. The polarisation of the e^- beams contributes differently to 
the various terms entering the total cross section, which is in the denominator of A^, and to the 
interference term ~ be in the numerator of A^. This results in a small net effect on A^ of the initial 
beam polarisation. This can also be inferred from the formulae of the various terms for polarised e^ 
beams which can be derived by using the formulae given in the Appendix. For higher cm. energies 
the asymmetry is larger, as the cross section in the denominator decreases faster with rising yfs as 
the numerator. However, the total number of events will decrease for the same reason, which can 
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Figure 9: Contour plot of the up-down asymmetry in the a — b parameter plane for Mjj = 120 
GeV and y/s = 800 GeV. The e ± beams are unpolarised. The grey code indicates the value of 



be balanced by a higher luminosity. 

The errors on a, b have been extracted in the same way as before with AO(a ,b ) of Eq. @ 
given by 



AO(a ,b ) = AA^ = -=Jl-A 2 . (21) 

The errors extracted solely from A^ turn out to be larger than the absolute values of a, b. This is 
the result of an interplay of a too small variation of A^ with a and/or b and too large errors AA^. 



We conclude this subsection by showing in Fig. [10] the up-down asymmetry in case we impose 
a 2 = 1 — b 2 . Being proportional to b it vanishes for 6 = 0. As it is also proportional to c, at b = 1, 
where a = 0, it vanishes due to our choice c = —a. As can be inferred from the figure, switching 
the polarisation from P e - = — 0.8, P e + = 0.6 to P e - = 0.8, P e + = —0.6 slightly increases A^. In 
this case, however, the cross section is not increased by about a factor of 2 any more as in the case 
of our original choice of P e ±, but is decreased by about 10%. 

2.5 The combined sensitivity 

As has been discussed in the previous subsections the total cross section is a good observable to 
determine a, whereas less good for the measurement of b. The top polarisation asymmetry, on 
the other hand, can only be used to distinguish a CP-odd from a CP-even Higgs boson. The up- 
down asymmetry A^ tests CP-mixing, but the errors on a and b are too large to be useful for a 
determination of a, b. The accuracy on a (b) will substantially improve, however, if the parameter b 
(a) has been extracted beforehand from some other measurement. We therefore combine all three 
observables Oi (i = 1,2,3) given by a, Pt and A^ to derive new sensitivity areas for a and b by 
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Figure 10: The absolute value of the up-down asymmetry as function of b with a? + b 2 = 1, 
for unpolarised beams (red/full), polarised beams with P e - = — 0.8, P e + = 0.6 (black/dash-dotted) 
and polarised beams with P e - = 0.8, P e + = —0.6 (green/dotted). 



performing a x 2 test, with x 2 defined as 

(Oil 

(AOi(ao,6o)) 



2 (OKM) -O^(ao,&o)) 2 , 00 s 
X = 2. fAO^ ^2 • ( 22 ) 



i=l,2,3 

The errors Aa 1 * 1 , Aft 1 * 1 are again determined by the maximal extensions of the insensitive areas and 



are shown at la confidence level in Fig. 11 for unpolarised and in Fig. 12 for polarised beams with 
y/s = 800 GeV and JC = 500 fb -1 . 

In case of unpolarised beams the plots essentially reproduce the errors on a and b extracted 
from the total cross section alone. This was expected since it is the most sensitive observable and 
Pt and do not help to improve the errors in this case. If the initial beams are polarised, however, 
the errors on a are remarkably reduced. This is due to the mutual interplay between a and Pt in 
constraining the parameter ranges of a and b. The polarisation is necessary in order to increase 
the cross section and Pt- For Pt this also leads to a smaller statistical fluctuation APt ~ lj\J aC. 
As for b, polarisation slightly helps to decrease the error for \b\ > 0.5. The observable A^ does not 
have any effect on the reduction of the errors at \/s = 800 GeV. 



At high cm. energies of y/s = 3 — 5 TeV as realized in the CLIC design 63 the total cross 
section, which scales with the energy, is smaller. On the other hand this will be balanced by higher 
integrated luminosities of J C = 3 — 5 ab _1 , so that the signal rate approximately remains the 
same. The top polarisation asymmetry does not change a lot compared to 800 GeV cm. energy, 
as can be extrapolated from Fig. [Tj which shows Pt as a function of \fs. The statistical fluctuation 
APt, which is inversely proportional to the signal rate, almost does not change either. The up-down 
asymmetry, however, gets larger with rising \/s. Therefore at high energies in the multi-TeV regime 
all three observables contribute significantly to x 2 - Their mutual interplay results in remarkably 
small errors on a less than about 0.2 in large parts of the parameter space. This can be inferred 



from Fig. 13, which shows for \fs = 3 TeV and C = 3 ab the errors on a and b at la confidence 
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Parameter a Parameter a 

Figure 11: Errors Aa + (upper left) and Aa~ (upper right) on a as well as Ab + (lower left) and Ab~ 
(lower right) on b, by combining all 3 observables a,Pt,A^, at la confidence level for M$ = 120 
GeV and y/s = 800 GeV with C = 500 fb _1 . The beams are unpolarised. The colour code 
indicates the magnitude of the respective error. 

level extracted from the combination of all three observables for polarised beams. It turns out that 
Pt and contribute to the reduction of the error in complementary areas of the a — b parameter 
plane. As for b, the error is < 0.3 except for low values of a and b values around ~ ±0.2, where it 
gets worse. 

3 Radiation of a spin 1 particle 

In this section we investigate the question to what extent the spin of the particle produced in 
association with a tt pair affects the total cross section and the polarisation asymmetry. The 
motivation is two-fold. On the one hand tiZ production contributes a large part to the irreducible 
SM background of associated Higgs production, so that a distinction of these two processes for 
the identification of the signal process is necessary |53j. On the other hand numerous models 
of New Physics predict the existence of additional neutral particles Z' with spin J = 1. In the 
following the associated production of such a particle will be investigated without referring to a 
specific modelj^ This particle is only demanded to couple to top quarks but not to e^, like the 

3 See Ref. |64| for associated SM Z boson production with anomalous couplings. 
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Figure 12: Same as Fig. 11, but for polarised beams. 



SM Higgs bosonj^] The description of the tiZ' coupling in an effective model-independent ansatz 
demanding Lorentz invariance and hermiticity and including terms of up to dimension 5 depends 
on ten parameters [65]. In a first approach we do not take into account anomalous vector couplings 
to ti of dimension higher than 4. This will be left to a future publication. The tiZ' coupling is 
hence given by 



a 



ttZ' 



ten/nig 



(23) 



with the vector and axial- vector couplings gy and qa chosen to be of 0(1). 

Assuming a particle with mass of 120 GeV has been produced at an e + e~ collider in association 
with a top-quark pair, it is investigated if the cross section and its energy dependence as well as Pt 
allow for the distinction of a SM Higgs H (and hence a spin particle) from a spin 1 particle. We 
therefore consider e + e~ — > tiZ' production of a Z' with mass Mz' = Mjjo = 120 GeV. The cross 
section is assumed to have the same magnitude as expected for the production of a SM Higgs with 
this mass value and a cm. energy of 800 GeV, within an error of 10%. This is the error on SM 
Higgs production to be expected for these parameter values at an ILC. This condition leads to a 
range V of values gy and gA, 



9V,9A G V 



vtiz'(gv,gA) = °tiH° ± io% 



(24) 



4 A non-vanishing coupling eeZ' would have been excluded by the LEP experiments for the masses we study. 
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Figure 13: Errors Aa + (upper left) and Aa (upper right) on a as well as Ab + (lower left) and Ab 
(lower right) on b, by combining all 3 observables a,Pt,A^, at la confidence level for M$ = 120 
GeV and y/s = 3 TeV with C = 3 ab" 1 . The e ± beams are polarised. The colour code indicates 
the magnitude of the respective error. 



For these values, the top polarisation asymmetry is calculated and compared to Pt(ttH°) = 0.11. 



Fig. 14 (upper left) shows the values of Pt(ttZ') for the range of values V . The white strips in the 
coloured area appear where Pt(tiZ') differs from Pt(ttH°) by less than 5a. These regions are small. 
We conclude that for a particle produced in association with a top quark pair with the same rate 
within the experimental error as expected from ttH° production, the top polarisation asymmetry 
can be exploited to distinguish a spin 1 state from a spin state if neither gy nor gA is close to 0. 

Moreover, the ratio of cross sections at different cm. energies can be exploited. This is shown 



in Figs. 14 upper right, lower left and lower right. The ratio of the ttZ' cross section at two different 
collider energies for three cm. energy combinations, respectively, is plotted for the parameter range 
V as defined above. All three ratio combinations are larger than the corresponding SM values and 
differ from them by more than 5a. The SM values are 

(W1000 GeV) <W1300 GeV) = 

a H ou (800 GeV) ' <W800 GeV) 

<W1300 GeV) _ 

WlOOOGeV) " °' 72 - (25) 
The SM ratios are very different from the ttZ' production ratios since the maximum of the SM 
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Figure 14: The Pt contours for ttZ' production with M%t = 120 GeV and yfs = 800 GeV (upper 
left). The contours for the ratio of the cross sections at different cm. energies are shown in the 
upper right plot for a(1000 GeV)/o-(800 GeV), the lower left plot for ct(1300 GeV)/cr(800 GeV) 
and the lower right plot for <r(1300 GeV)/cr(1000 GeV). The colour code indicates the magnitude 
of Pt and of the ratios, respectively. The couplings gA and gy have been chosen as described in the 
text. 



cross section occurs at ~ 750 GeV and hence all considered ratios in case of the SM Higgs are less 
than 1. For tiZ' production the maximum lies at 1000 GeV or higher, depending on the values of 
gv, gA- Although these results of course depend on the mass of the particle radiated from the top 
quark pair and on the considered energies, they show that the top polarisation asymmetry and/or 
ratios of cross sections are sensitive to the spin value of the particle produced in association with 
tt. 



This statement is confirmed by the plots presented in Figs. 15 and 16 Figure 15 shows the 



total cross section a(tiZ') compared to the cross section of scalar Higgs production in association 
with a top quark pair, a(ttH). The H and Z' masses have been chosen equal, Mz> = Mh = 120 
GeV. For the Higgs coupling to the top quark pair the maximum coupling strength a = 1 has been 
adopted. The couplings gv and gA have been chosen within a parameter range V' such that the 
maxima of the two cross sections in comparison are the same, 

gv,g A eV' <=> a^(g v ,g A ) = a^ . (26) 

Three cases are studied, where either gy or gA are non-zero and where both gy and gA are non-zero. 
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Figure 16 shows the comparison of a(tiZ') with the associated production of a pseudoscalar A 
with a top quark pair, a(tiA), for Mz 1 = M A = 120 GeV. For the CP-odd coupling parameter 
to the top quarks 6 = 1 has been assumed. The coupling parameters gy, gA have been chosen as 
before within a parameter range V" such that the maximum values of the cross sections are the 
same, 

9v,9AeP" a^(g v , 9A ) = a^, (27) 

and once again three types of coupling combinations gy and gA are investigated. As can be inferred 
from the figures, the curves for spin 1 production and CP-odd Higgs production are very similar 
both in the threshold rise and in the continuum. This is to be contrasted with the comparison of 
the curves for spin 1 production and CP-even Higgs production. The threshold rise is not so much 
different. But in the continuum they clearly differ. This reproduces the difference in the ratios of 
cross sections observed above. 



4 Summary 

Once a Higgs boson has been discovered at the LHC, its CP properties have to be determined. 
Should the Higgs boson not be CP-even but CP-odd or a CP-mixed state, then the top sector 
represents an ideal working ground to test its CP properties, since the various possible CP-parts 
of the Higgs boson couple democratically to the top quarks. This is in contrast to the Higgs 
coupling to gauge bosons which projects out the CP-even part of the coupling and has only a small 
admixture of CP-odd parts from loop diagrams. In this paper, we investigated the production of 
a spin state with arbitrary model-independent CP properties in association with a top quark 
pair at a future e + e~ linear collider. The CP properties of the Higgs coupling to the top quarks 
have been parametrized in a model-independent way by a parameter a for a CP-even Higgs, by a 
parameter b for a CP-odd Higgs and by simultaneously non-vanishing a and b for a CP-mixed state. 
These parameters can be determined by a measurement of the total cross section, the polarisation 
asymmetry of the top quark and the up-down asymmetry of the antitop quark with respect to the 
top-electron plane. The former two observables are CP-even and can be exploited to distinguish 
a CP-even from a CP-odd Higgs boson. Since the up-down asymmetry A^ is CP-odd, it can be 
exploited directly and unambiguously to test CP violation. 

The sensitivities to a and b have been studied in each observable separately before investigating 
the combination of all three observables. We found that the total cross section is most sensitive to 
a and to some extent to b. The observables Pt and A^ do not exhibit much sensitivity to a and b, 
although polarisation of the initial beams slightly improves the sensitivity in case of Pt- The 
combination of all three observables, however, remarkably reduces the error on a for polarised 
beams. At 800 GeV cm. energy this is due to the mutual interplay of a and P% in constraining the 
parameter range of a and b. At multi-TeV energies A^ is larger, and it is now the interplay of all 
three observables which further reduces the error. 

We also investigated to what extent the radiation of an arbitrary spin 1 particle from the top 
quark pair, which like the Higgs boson does not couple to , can be distinguished from a spin 
state with the same mass. It turns out that the top polarisation asymmetry and ratios of cross 
sections at different collider energies represent good observables to tell a spin from a spin 1 state. 

In summary, associated Higgs production with a tt pair at a future e + e~ collider allows for 
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g v = 1.8, g A = 0, a = -c = 1, b = 






Figure 15: Total cross section a(tiH) (red/full) compared to a(ttZ') (green/dashed) as function 
of the cm. energy for Mz> = Mh = 120 GeV. The parameters gv / (upper), gA ^ (middle) 
and gv,9A 7^ (lower) have been chosen such that the maximum values of the cross sections are 
the same. 
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Figure 16: Total cross section a(ttA) (red/full) compared to a(tiZ') (green/dashed) as function 
of the cm. energy for Mz> = Ma = 120 GeV. The parameters gy ^ (upper), gA ^= (middle) 
and gv,9A 7^ (lower) have been chosen such that the maximum values of the cross sections are 
the same. 
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the determination of the spin and CP properties of the Higgs boson in a model-independent way. 
Combining the total cross section and its energy dependence, the top polarisation asymmetry and 
the up-down asymmetry considerably helps to reduce the errors on the coupling parameters of the 
Higgs Yukawa coupling which reveal the CP nature of the Higgs boson. 



Appendix 

Cross sections for initial beam polarisation 

The cross sections oq and o\ as defined in Eq. ^ get modified for initial beam polarisation. The 
polarisation acts differently on the contributions to the cross sections from the various diagrams. 
The changes can be summarized by replacing the coupling factors in the cross sections. Denoting 
by P e - and P e + the degree of longitudinal polarisation of the electron and positron, respectively, 
the changes to be applied for the couplings appearing in a = 2ao Eq. Q and <7i Eq. (12) are 

QlQt -»• Q 2 e Ql{l-P e+ Pe-) 

2QeQtv e -> 2Q e Qt[v e {l-P e+ P e -) + a e {P e+ -P e -)} (28) 
{v 2 e +aD ->• [(v 2 + a 2 e )(l-P e+ P e -) + 2v e a e (P e+ -P e -)] . 

The up-down asymmetry originates from the interference term between the diagrams where the 
Higgs is radiated from the top line and the one where the Higgs is radiated from the Z boson. The 



denominator of As is simply given by the total cross section cf. Eq. (20). The differential cross 



section for the interference term ~ be reads in terms of the scaled Mandelstam variables S\ , S2 , T\ , T2 

dabc 3a 2 be gttH9zzH\fh t / s 

dSxdSzdT^ ~ 8tt 2 s3 h z (Si - h z ){S 2 - h){-l + S 1 + S 2 - h - h*) X 

r QeQtatae n $ 2a e v e a t v t # a 2 t {a 2 e + v 2 e ) * •> PlP2<laqb , . 

i Jh^i) Dl + Th^W 2 + (h z - 1) 2 3 / ' (29) 

with hx = fn 2 x I s an d the form factors given by 

Df = (h* + 1 - Si)(h* - 2h z - 1 + Si) (30) 
D% = (l-Srf-hl (31) 
Df = (-/i$-2/i t -l + 5i + 2 < S 2 )(/i$-l + 5i-2r 2 ) . (32) 

The scaled Mandelstam variables in terms of the reduced energies x\,X2, the azimuthal angle x 
of the antitop, the polar angle 9 of the top quark, the angle 9\ 2 between top and antitop and the 
velocities fii } 2 = yfl — Aht/x 2 2 read 

S\ = 1 — /i$ — X3 (33) 
S 2 = 1 + ht-xt (34) 

Ti = ht- ^i(l -A cos 0) (35) 
T 2 = /i$ — -xi[x3 — sin 6 sin 9\2 cos x x %fii ~ cos 6( x iPi + cos 612)] ■ (36) 
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In case of polarised beams, the following replacements for the coupling factors have to be made 

QeQta e -> QeQt Ml - P e +P e ~) + V e (P e+ - P e -)] 

2a e v e -> [2v e a e {l-P e+ P e -) + (a 2 e + v 2 e ){P e+ - P e -)] (37) 

(«e + «e) -> [(«e + t£)(l " ^e+^e-) + 2a e ^ e (P e+ - P e -)] . 
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